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Abstract 

Starting with the two-loop effective potential of the MSSM, and assuming a super- 
symmetric scale well above Mz, we derive a simple analytical approximation for the 
lightest CP-even Higgs boson mass including resummation of higher order logarithmic 
terms via RG-improvement and finite non-logarithmic terms up to 0(a s a t ). This for- 
mula describes the most relevant radiative corrections to the MSSM Higgs boson mass, 
in particular, those associated with non-zero top-squark mixing. 
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1 Introduction 



Supersymmetry (SUSY) stabilizes the hierarchy between the scale of electroweak symmetry break- 
ing and the fundamental energy scale (the GUT, Planck or string scale) and makes plausible that 
this breaking occurs in a weakly coupled Higgs sector. The simplest realistic model that accom- 
modates these ideas is the Minimal Supersymmetric Standard Model (MSSM), the paradigmatic 
testing ground of low-energy SUSY. Given a weakly coupled Higgs sector and a scale of electroweak 
symmetry breaking v = 246 GeV, it follows [p]] that the spectrum of the theory must contain an 
scalar particle h° with mass controlled by the Fermi scale and non-zero couplings to the W and 
Z° gauge bosons (which is crucial for its detection at accelerators ||). In the MSSM, the mass 
of this light Higgs boson (besides this Higgs particle, the Higgs sector in the MSSM contains a 
heavier scalar H°, one pseudoscalar A and a pair of charged Higgses H ) is calculable. A precise 
determination of this mass is of paramount importance for the experimental search of SUSY. If the 
experimental lower bound on the Higgs boson mass || increases above the theoretical prediction 
we can rule out the MSSM. It is therefore understandable that this calculation has received a great 
deal of attention and attracted the efforts of many groups j4j|-[fL0|. 



At tree level, the mass squared, m\o, of the light Higgs boson has an upper bound (saturated 
for large values of the pseudoscalar mass m^o) given by Mf cos 2 2j3. This is already below the 
experimental lower bound from LEP2 ||. However, radiative corrections can raise the upper 
bound on m\ dramatically ||. The dominant contribution is 



3mf , m 2 



Am 2 h0 = -^ln— | , (1) 
2n 2 v 2 mi 



where m t is the top quark mass and mi a common top-squark mass. Radiative corrections to m 2 h0 
have been computed using different techniques: effective potential approach, direct diagrammatic 
calculation and effective theories with renormalization group (RG) tools. Each approach has its 
own virtues. The effective potential way simplifies the computation; the diagrammatic calculation 
is unavoidable to pick up some particular corrections; and the RG approach permits resummation 
of logarithmic terms to all loops and provides a physically meaningful organizing principle for the 
radiative corrections. There is no real need to choose one among these methods: the best way to 
proceed is to combine the three methods, taking advantage of the best virtue of each in turn. 
The n-loop contribution to the (dimensionless) quantity m 2 h0 /m 2 has the schematic form 

±(aln^) k a"-\ (2) 



fc=0 



m t J 



where a represents the expansion parameter (with at = h 2 /4ir and a s = g 2 /47r giving the dominant 
contributions), and the logarithms \n(mi/mt) can be sizable for large values of m^. The term k = n 
in Eq. @ is the leading-logarithmic correction and dominates the n-loop contribution for large 
m~ t . The term k = n — 1 is the sub-leading logarithmic term, etc. Finally, the term k = 



1 



gives the finite non-logarithmic piece of the n-loop correction. The complete one-loop radiative 
corrections to m 2 h0 /m 2 have been calculated. The dominant leading- logarithmic contribution is 
the atln(raf/rat) term in Eq. ([!]). At this loop order there are no a s corrections, which enter 
at two-loops. The most important finite corrections depend on the top-squark mixing parameter 
X t = A t + fM cot f3, is given by 

. 2 3mj (X 2 X? \ 

and its effect on m\o can be important. The maximum value for the upper bound on the Higgs 
mass is obtained for X t 2 = 6m? (the so-called 'maximal- mixing' case). 

The most important part of the higher order radiative corrections can be collected and re- 
summed using renormalization group techniques f|. Resumming with one- loop RG equations 
takes into account the leading-logarithmic corrections to all loops, and using two-loop RG equa- 
tions sub-leading logarithms are also included. However, the finite non-logarithmic terms cannot 
be obtained in this way. Getting them at two-loops requires a genuine two-loop calculation. This 
was first done, in some particularly simple limit, in Ref. 0, but the effect of non-zero top-squark 
mixing was not included. Recently, a two- loop diagrammatic computation to 0(at<y s ) has been 
completed ||, and the effective potential has been also calculated to the same order 0. These 
studies show that the most significant two- loop effect, not previously taken into account by RG 
techniques, comes from the finite pieces dependent on the top-squark mixing. The new refined 
upper bound on m^o in the region of maximal mixing can increase up to 5 GeV (if m 4 ~ ~ 1 TeV) 
and the condition for maximal-mixing itself gets also slightly modified.^ 

Making good use of RG resummation, it is possible to obtain compact analytical approxima- 
tions for m\ which take into account the most important radiative corrections ||. Besides being 
of direct practical interest, these formulae are theoretically interesting as they provide a clear 
picture of the physical origin of the dominant contributions. In this paper we extract such an 
analytical approximation for m\ starting with the two- loop effective potential computed in f9|. 
Our final formula includes the most important 0(a s a t ) radiative corrections, in particular the 
finite terms associated with non-zero top-squark mixing. Our results agree with those obtained 
by diagrammatic techniques |T(| where they overlap, but are computed by an alternative way. We 



include, in addition, RG resummation which allows us to write a particularly simple final formula 
for the radiative corrections to m 2 h0 \ this is the main result of our paper. Again, to the order 
at which previous RG results || were computed, we find agreement with our results. Our final 
formula for m 2 h0 improves over previous RG formulae by including the genuine two-loop threshold 
corrections (which are important for large values of the top-squark mixing) and over previous 
diagrammatic results by incorporating RG-resummation of logarithmic corrections. 

1 This is not a disagreement between different calculations, but the result of including in the analysis of M 
two-loop effects inaccessible to previous approaches. 



2 



2 Two-loop effective potential in the MSSM 



The two-loop effective potential in the MSSM to 0{a s ott) has been calculated for the case of zero 
top-squark mixing and tan/? = oo in 0, where tan/5 = v 2 /vi is the usual ratio of the vacuum 
expectation values of the two Higgs fields; and for arbitrary top-squark mixings and tan (3 values 



m 



(16tt 2 ) 2 V 2 {h x , h 2 ) = Sg 2 3 U(m t ,m t ) -2m, 2 l(m t ,m t ,0) 

\ 2 2 

A i=l i=l 

2 r 

+2J(m § , m t ) - J2 J( m u, m s ) + J(m~ u , m t ) + (m| - m? - m 2 t )I{m~ u ,m~ g , m t ) 
i=i L 

-4m § m t StQ J(m fi , m § , m t ) - I(m~ t2 , m § , m t ) 



(4) 



where the arguments hi, h 2 are the neutral scalar components of the Higgs fields Hi and H 2 . All 
the masses and mixing angles in Eq. (f|) should be understood as hi, /^-dependent quantities, for 
example, the top quark mass m t = -^h t h 2 . We use short-hand notations c t = cos^, s t = sin6^, 
where 9i is the top-squark mixing angle. The minimally subtracted^ two-loop scalar functions / 
and J are 



;(m 1 ,m 2 ,m 3 ) 



2 2 

/ o 2 2\ i ^^2 i ^7^3 

v — m l + m 2 + m 3 ) m-^r In 



q 2 q 2 
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+ (mj - m 2 + m 3 ) In — In — + (m x + m 2 - m 3 ) In — In — 

m 2 2 2 

- 4(m 2 In — | + m\ In — | + ml In — |) + ((mi, m 2 , m 3 ) + 5(m 2 + m\ + ml) 

q 2 <y 2 <y 2 



J{mi,m 2 ) = m\m\ 



1 2 2 2 

, mf , m- , mf , m, 
ln ^- ln # + ln # ln # 



(5) 
(6) 



where Q is the renormalization scale. The function ((mi, m 2 , m 3 ) in the expression of J(m 1? m 2 , m.3) 
has been calculated in fl2| by a differential equation method and in [H| by a Mellin-Barnes inte- 
gral representation method. Its final form can be expressed in terms of Lobachevsky's functions 
or Clausen's integral functions and their analytical continuations. 



2 We choose to work in the DR scheme jTT|, as required by preserving the SUSY Ward identities. This choice 
also simplifies complications associated with vector-boson loops in the MS scheme. The two-loop effective potential 
in Eq. (Q) is calculated in the Laudau gauge, but it is also correct for a general gauge. 

3 Subtraction of one- loop sub-divergences can be treated following the first reference of In obtaining Eq. (Q), 
we have also added freely one-loop diagrams depending only on m§, as they do not affect the Higgs boson mass 
calculation in the effective potential approach. 
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It is important to check that the effective potential V(hi,h 2 ) = Vo(/ii,/i2) + ^1(^1,^2) + 
V 2 (hi, h 2 ) is invariant under changes of the renormalization scale up to higher order terms of 0(a 2 ). 
To prove this, we write the tree-level and one-loop potential as follows (neglecting contributions 
from the gauge boson, Higgs/Goldstone boson and neutralino/chargino sectors): 



Vb+i(/iiA) = \{^HA^) h l + \( m H 2 +^) h l+ B ^ h i h 2 + ^ G{m h )+G(m k )-2G{m t ) 



, (7) 



where m Hl , m H2 and £> M are the soft-breaking Higgs sector mass parameters, /i the supersymmetric 
Higgs-boson mass parameter, and 

G(m) = —(in—--) . (8) 
It is then straightforward to show that 



2 V Q'' 



modulo terms independent of the Higgs field h 2 , and where X t = A t + /icot/?, with A t a tri- 
linear coupling (with dimensions of mass) appearing in the following term of the soft-breaking 
Lagrangian: h t A t H 2 QU . and are soft-breaking masses for the left- and right-handed top 
squarks, Q and U, respectively, and M 3 is the gluino soft mass. [It is a nontrivial check that all 
InQ 2 terms cancel with each other in Eq. (B|)]. In the above equation T>^ and T>^ stand for one- 
and two-loop RG variations of the parameters, and we have used the following equations 

dI{m ^ m) = m ;in^ + m |.n=| + m |.ng-2K +m I + m 3), (10) 

dlnQ 2 = m ^{ 2 - ln Qt- ln W ' 
and the MSSM RG equations from Ref. p|, 



dm% 3h 2 



.m 



i^+Ml + Ml + A 2 ) 



dlnQ 2 16tt 2 

16 * 2 (m 2 H2 + M 2 ~ + Mf f + A 2 + 2M 2 - 2M 3 A t ) , (12) 



(16tt 2 ) 2V m Q u 
<91n/i dlnh 2 3h 2 t 8gjh 2 t 



dlnQ 2 dlnQ 2 32vr 2 (16 

dB u 3h 2 (Bu t \ ISglh 2 fB„ , , , 



ainQ 2 16tt 2 V 2 "7 (16vr 2 ) 2 V 2 

We will see that the RG does not only provide an important check for the Higgs boson mass 
correction formulae in the later sections, but also allows us to present those formulae in a more 
physically appealing form— the RG improved form. 
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Using the effective potential in Eq. (|), two- loop radiative corrections to m 2 h0 [to the order of 
0(a s a t )} have been calculated numerically in |§. In the following, we shall derive an approximation 
formula valid for M^{= M^) and m^o 3> Mz, where m^o is the mass of the pseudoscalar Higgs 
boson A . 



3 Analytical expression for Amjo from the effective poten- 
tial 



To simplify our analytical expression, we assume the squark soft masses satisfy Mq = = Ms, 
where Ms is the SUSY scale. The two eigenvalues and mixing angle of top-squark squared-mass 
matrix are 



m,- 



m~ + m t X t 



m 



t2 



mf — m t X, 



t , 



1 

72 



(15) 



where the average top-squark mass m~ = Mj + m 2 . We shall further assume the gluino mass 
nig = M s - 

Under these simplifications, the effective potential to 0(a s a t ) two-loop order is (neglecting 
one-loop sub-dominant terms) 



V{h u h 2 ) = V {ht, h 2 ) + 



16tt 2 



2G(m { ) - 2G(m t ) + m 2 X 2 In p ™ W 



12m 4 



+ 



%1 



(16vr 



2^2 



J(m t , m t ) — 2mfl(m t , m t , 0) + J(m 4 -, m^) + 2m ( -/(m ( -, m~ t , 0) 
+2 J(M S , m t ) - 2J(M Sl mi ) - 2 J(m t , m f ) 



+ ™t 



mi 



2X t m,(l-ln^|) 2 + X, 2 (l-ln 2 ^|. 



mi 



m 



mi 



-± (l-ln-M 2 + 21n^ + + 



Q 



mf 



X? 



mi 



3m? 



mi 



21n gt 



12m 4 



•,(16) 



where we have separated terms depending on X t (to make the top-squark mixing effects transpar- 
ent) and expanded the effective potential in powers of mt/mi (and mtX t /mf) with higher order 
terms in m t jm- t neglected. It is now straightforward to find the corrections to the four entries in 
the Higgs boson squared-mass matrix 4 



AM 2 n 



d ( 1 dV{h u h 2 ) 



dhi \hi dhi 



/ hi=vi,h,2=V2 

4 The effective potential approach evaluates these corrections at the zero external momentum limit. This is 
consistent with our other approximations since the tree-level light Higgs boson mass is related to the electroweak 
gauge couplings and can be neglected. 
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AM 2 12 = 



AM 



22 




2 

m~ 



-flA t 111 g| + 



m\[iA t X\ 
6m~ 



m x 



m~ t \ 1 



mill 
+ — ^tan/? 
2m- t 

, rn*X t 
m| 



m 



9 

(1 - In— f) 2 + 21n 



2\ 2 



A 1 



1 2 m ? 
Q 2 



fiXt 



Q 2 

tan/?) 



mf 
ml 



m 2 nA t 



m~ 



tan/? 



2k # 



+ 



2m| 



-/i H tan/? 

3 



AM 2 21 = 



d 2 V{h u h 2 



dh 1 dh 2 



hl=Vl,fl2=V2 



3oj 

47T 



m 



m*iiA t X^ 



mf 



2mf 



m|/i 
m t ~ 

,2 - 



mA 1 — In 



m 2 x 2 



Q 2 



+ A t 1-ln 2 



1- 51n — ^ n 2 -4 

2 g 2 2 g 2 



mf 

g 2 

2 n 



31n— | 

mf 



+ 



mffi 



mf 



A t + 2X t 



miX; . , uX , 



mf 



i- 2 - 

, mf 


^Wt 2 1 


. 1 " 21n #. 


2mf / 



/l 2 



9 / 1 ov(h u h 2 y 



dho \ h 2 



dh 2 



hl=Vl,h2=V2 



^^m 2 In — % + — < 3 cot 8 
- 4ir{ 



m 2 



2a s a t 2 
H —m t 



7T Z 



i 2 m ? 

In 2 -4 

mj 



. , mf mfiiA t X. 

, 9 m~ , 9 m? , m~ 
2 In 2 -1 + 2 In 2 + In -4 - 



+ 



12 AX, 



mf 



A 2 ^ 
mj 



—fi cot /? 



m.- 



mi 
+— 

mi 

,2 



^cot/?((l-ln^|) 2 + 21n^|) +2(A + X t )(41n^|-21n^|) -4A 



Q 2 

mf 



g 2 



i 



"it 



+ A(l-ln 2 ^|) 



g 2 



mf 



m 2 A 



mf 



Q 2 



AX t - fx cot - 8X t In 



m 2 



mf 



mf 



mf 

g 2 



+ 



mjX t 
mf 



g 2 



^-^cot/? + ^ 



m 2 



(17) 



(18) 



1 -2 In 



m 2 f 



Q 2 



m 2 A t X 2 
2m\ 



X t 



(19) 



where a s = g 2 /47r and a t = h\j^. Again we have neglected terms in higher orders of m t jm- t . 

The Higgs boson mass to the two-loop order can be obtained by diagonalizing the CP-even 
Higgs boson squared-mass matrix M 2 by including these corrections 



M' 



Mf cos 2 8 + m 2 A0 sin 2 8 + AM 2 n -(Mf + m 2 AQ ) sin 8 cos 8 + AM 2 12 
-(Mf + m 2 A0 ) sin (3 cos (3 + AM 2 21 Mf sin 2 (3 + m 2 A0 cos 2 8 + AM\ 2 



(20) 
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This means that the mixing angles at two-loop order are different from those at tree-level. Nev- 
ertheless, this difference is small to a good approximation and we can therefore use the tree-level 
mixing angle a to compute the Higgs boson mass. We further simplify the mass correction formula 
by approximating cos a = sin/? and sin a = — cos/?, which are valid in the limit m^o ^> Mz- The 
final analytical expression for the two-loop order correction to the lightest CP-even Higgs boson 



mass is 



Am 2 



3m; 



h° 



2n 2 v 2 

+ — 

mz 



In— ^ + -4 

mf m~ 



1 + 2to # 



Ylm\ 

2 

, m, 

21n # 



+ 



2 

m x 



m. 



m 2 



tt 3 v 2 



\n z -4 - 2 In 2 — f - + 2 In 2 — £ + In -4 - 1 



+ 



^ 2 



3m? 



Q 
i -i 



Q 2 
xt 

12m\ 



nil 



(21) 



We note that all parameters in Eq. (|2T| ) are running parameters evaluated in the DR -scheme and 

ioi 

up to terms of order 0(a 2 ) by using the following equations 



satisfy MSSM RG equations. It is also easy to check that Am^ is indeed two-loop RG invariant 



d In m? 



din m 2 



d\nQ 2 «91nQ 2 



4a, 
' 3tt 



91nQ 2 



4a £ 



M 3 



(22) 



with the gluino soft mass M3 = = Ms in our approximation. 

At first sight, Eq. ([H]) seems quite different from those formulae obtained by the RG improved 
one-loop effective potential approach [Q and by the two- loop diagrammatic approach JTIJ. The 
difference can be settled by observing that the parameters in Eq. fl2~lD are MSSM running parame- 



ters, while the top quark mass in || is a SM running mass and X t and m- t in Jl0| are on-shell (OS) 



parameters. Our one-loop parameters receive radiative corrections from the SM as well as SUSY 
particles, so they are different from the parameters in other approaches where the SUSY particles 
are explicitly decoupled. Furthermore, part of the difference arises from converting parameters in 
the DR -scheme to the MS -scheme. 

To demonstrate this point, we first compare our result in Eq. ( pT|) with that of Ref. ||10|| . We 
use the following relations of DR top quark mass (mt) and top-squark running mass (mz) (in the 



MSSM) and the MS running top quark mass {m t ) and OS top-squark mass 
which for example can be deduced from Refs. |7|, [L4| 

m t (Q) = m t (Q) 

m { (Q) 



mz 



3^ 



m~ 

1 + ln # 



Xt 
mz 



mz 



4a. 2a, 



+ 



37T 37T 



2 

mi 
In— *■ 
Q 2 



[in the SM), 

(23) 
(24) 



up to higher order terms in mt/mz!^ All terms in Eqs. (|23D and 
top squarks and gluinos, except for the finite term — a s /3ir in Eq. (|2 



come from the decoupling 
which is a DR -MS -scheme 



^Parameters inside the brackets can be running or OS ones, as the difference is of higher order. 
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conversion factor. Similarly, one can relate X t between different approaches by considering the 
top-squark/gluino contributions to the self energy Ilf^ 0, which gives 



x t {Q) 



x t + 



m t 



KeU~ tL - tR {p 2 



X t ReZ t (p 2 



m. 



3vr 



mi 



m~\ ( m~\ X?' 

4 In -I: \+X t 5 + 31n-| - -L 



Q 2 



ml 



mi 



(25) 



where X t is the OS top-squark mixing parameter. The external momentum of ^-i L t R nas been set 
to the average top-squark mass TOfQ The top-quark self energy S 4 is evaluated in the MSSM and 
contains QCD corrections from top-squark/gluino and top/gluon loops (it can be deduced from 
Refs. 
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Re£ t (p 2 



m. 



a s 
Sir 



,m 2 , m~ X t 
31n-i + ln-f --^ 



Q 2 



Q 2 



m; 



(26) 



We note that in the effective SM the average top-squark mass m t - and the parameter X t are frozen 
and equivalent to their physical 'pole' values since all SUSY particles have been decoupled. 
Substituting Eqs. (0), (0) and (|25|) into Eq. pi]) , we obtain the mass correction formula 



Am 2 



h° 



3mf 
2n 2 v 2 

6X t 
+ — - 

mi 



\ ml X 2 1 Xt 

In ■ = r + ^ 

m\ m~ 12 m~ 



a s mf 
ir 3 v 2 



2 2 2 2 

-4 - 3 In 2 ^ + 3 In 2 ^| + 3 In - In ^ 



+ 



XI 



m 2 



2 r 



2 2 

, mf mj 

31n #- 31n # 



+ 



4 r 



Q 2 

1 , m?" 1 
12 + 2 ln ^ 



mf 
Q 2 



17 



(27) 



where we have explicitly indicated the definition used for different parameters in the one-loop cor- 
rections; different definitions for the parameters in the two-loop corrections would give differences 
of higher order. This equation is identical to the mass correction formula in [nj for Q = m t . 
Comparison of our formula with the one obtained using the RG-improved one-loop potential f| 
is left for the end of next section. 



4 RG-improved Higgs boson mass 

The final expression for m 2 in the last section was very convenient to make contact with the 
diagrammatic results of |l(J. However, a simpler and more transparent expression can be obtained 
by transforming to the RG language. Furthermore, the improvement of the formula goes beyond 
purely aesthetic reasons, as it resums higher order corrections as well. 

The idea is to let all parameters in the formula for m 2 h0 be running parameters and to choose 
the scale at which they are evaluated in such a way that higher order logarithmic corrections are 
automatically taken care of. Moreover, the scale at which each parameter has to be evaluated to 
achieve this, is susceptible of physical interpretation. 

6 Choosing p 2 = m~ ± m t X t changes the result by terms of higher order in m t /m^. 
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Starting from Eq. (pi]), we first use Eq. (^) to translate the MSSM top-quark running mass 
m t (Q) into that of the SM, m t (Q). We then use the solutions to the SM top quark mass RG 
equation and the MSSM RG equations of mi and X t in Eq. (|22|), 



m 2 t (Q) 


= m 2 (Q') 


L 2a Q' 2 ' 

7T Q 2 _ 


(28) 


mf(Q) 


= mt(Q>) 


i 1 + 3vr ln q 2 


(29) 


X t (Q) 


= X t (Q')- 


3vr 3 Q 2 


(30) 



to relate their values at two different renormalization scales, Q and Q'. Our final formula is 



Am 2 



h° 



2tt 2 v 2 



rr4(Q t )ln 
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thJ 



m?(Q th ) UrnHQ 



th) 



+ 



a s mf 


2X 4 


* t 2 


7X ( 3 

h + 


1 xt 


ix t 5 


7T 3 f 2 


m t ~ 


mf 


3 m? 


12 


6 



(31) 



th 



with Q\ = m t mi, Q\m t = Q'\ (satisfied by Q~ t = Q' t = m t or Qi = m t ~, Q'^ = m t mf) and Q 
Note how nicely everything falls into place. All higher order logarithms are reabsorbed into the 
one-loop RG-improved term and the only two-loop pieces remaining are the finite X r dependent 
terms. 

The logarithmic term in Eq. ([H]) can be interpreted as the result of integrating the RG equation 
for the quartic Higgs coupling between the high energy scale mi and the top-quark mass scale m t 
in the SM, which is the effective theory below mi. The coupling in front of this logarithmic term 
has to be evaluated at the intermediate scale Qt = ^/m t mi, choice which automatically takes into 
account higher order effects. Having computed the two-loop results, we can also say something 
on the scale at which the masses entering the logarithm should be evaluated, although not in an 
unambiguous way. If we insist in evaluating them at the same scale, that scale turns out to be 



Qi = Q't = m t- If; on the other hand, we prefer to keep Qi 



mi 



then Q' t = (m t m 2 ) 1 ^. The 



finite non- logarithmic 'one-loop' term is interpreted as a threshold correction for the quartic Higgs 
coupling at the scale mi, at which the SM is matched to the MSSM. In accordance to this, all 
the parameters in this term are evaluated at the threshold scale Q t h = mf- Then, the remaining 
two-loop finite terms give a two-loop contribution to this threshold correction.^] Eq. ( |5T| ) can be 
slightly modified to automatically include the logarithmic terms of order 0(af) (and higher) if we 
refine Eqs. (p3|), (p8|), (|29"D and (|30|) to take into account 0(a t ) corrections. 



7 Actually, an alternative way of computing this two-loop threshold correction is the following. Start with the 
expression for the effective potential in Eq. (f 6) and make use of Eq. ( p3| ) to trade between the running masses in 
MSSM and SM. Expand then the potential in powers of mt- In that expansion, the finite term of order m\ gives 
the correction to the Higgs quartic coupling, in accordance with our final result Eq. 
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Figure 1: Higgs boson mass m h o versus X t /m^. The solid lines and dotdashes correspond to the 
two- loop Higgs boson masses calculated from the RG-improved mass formula Eq. (|3ll) , with and 
without including the two-loop finite threshold corrections. One-loop masses are shown in dashes 
for reference. 

We show in Fig. [I] the numerical result from the RG-improved Higgs boson mass formula Eq. 
(|3ll). The two- loop finite threshold corrections are generally small for small and moderate mixing 
parameter X t /rhi- For large mixing X t /rhi > 2, however, the two-loop threshold corrections can 
contribute corrections of the size of ~ 3 GeV. 

In terms of the parameter X t (Q t ^) fm^Qt^) , we can obtain the refined condition for maximal 
mixing as 



X t (Q 



thJ 



m- t {Q 



thJ 



±\/6 + 



5a s 



(32) 



which corresponds to a 2% shift to the position of maximal mixing. This condition can be rewritten 
in terms of OS quantities as 



(33) 



We are now ready to compare with the formulae presented in || using the RG-improved one- 
loop potential (some two-loop input was used in the last paper in ||). For zero squark- mixing, 





= ±V6 


a, / , m?\ 
1 - — 3 + In -| 


a s 
+ — 


m~ t _ 


max 


7T V mi) _ 


7T 
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the leading and next-to-leading logarithmic terms in Eq. ( p7|) agree with || but those papers also 
gave the RG-improved one-loop threshold correction for non-zero top-squark mixing case, which 
is [omitting the 0(af) terms] 



3rnf (mt) 



2vrV 



X, 



m 



7T 



m| 

2 

mi 



(34) 



In this formula the X t parameter is implicitly evaluated in the MSSM at the threshold scale, 
exactly as in our formula Eq. (0), and agreement between the two is found after expressing 
m t (m t ) in terms of m t (Q t ^). However, the two-loop finite terms included in Eq. (|3T| ) can not be 
reproduced by the RG-improved one-loop effective potential approach. 

To summarize, our analytical expression for m\ agrees with those obtained by explicit two-loop 



calculation |10[ and RG-improved one-loop effective potential approach |j| where they overlap. We 
note that parameters in different models {i.e., MSSM and SM) and renormalization schemes have 
been used in previous studies, our discussions should have resolved any possible confusion. 



5 Conclusions and outlook 

We have used the knowledge of the MSSM effective potential up to two-loop 0(a s a t ) to extract a 
simple analytical approximation to the mass m h o of the lightest CP-even Higgs boson in the simple 
case of a single supersymmetric threshold significantly higher than Mz- We have considered the 
case of arbitrary tan/5 and non-zero mixing in the top-squark sector. We have derived a RG- 
improved formula to resum large logarithmic corrections, achieving a particularly simple and 
illuminating final result. 

Our results agree with previous analyses based on the RG-improved one-loop effective potential 
and diagrammatic calculation to the order at which such studies were performed. By doing 
this we clarify the relation between different approaches and identify the two-loop origin of the 
discrepancies between || and || for large values of the top-squark mixing. This difference can 
be attributed to a two-loop threshold correction, most easily calculable in the effective potential 
approach, as we have shown. We emphasize that in applying two-loop mass correction formulae, 
one should be particularly careful about the parameters entering the one-loop formulae; we have 
done this by differentiating running and OS parameters in our mass correction formulae, Eqs. 
(H and (0). 

The combined use of the three different techniques (effective potential, diagrammatic calcula- 
tion and RG-resummation) is very powerful and should be applied to compute the still missing 
0(af) corrections to m 2 h0 , which can be similar in magnitude to those analyzed in this paper. 
Phenomenological analyses done with the expressions currently given in the literature are accu- 
rate only up to the inclusion of these 0(af) corrections, which one could expect to give a shift in 
m h o not greater than 5 GeV, but quite interestingly, we expect that this shift goes in a direction 
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opposite to that of the two-loop 0(a s a t ) corrections. The leading and next-to-leading logarith- 
mic corrections may be obtained by a RG resummation approach, but the finite non- logarithmic 
corrections can only be extracted by a direct two-loop calculation to the order 0(af), e.g., from 
effective potential or from explicit diagrammatic calculation. 



Note Added 



A comparison between the diagrammatic and RG results for m? is mentioned in M. Carena 



S. Heinemeyer, C.E.M. Wagner and G. Weiglein, ||hep-ph / 99 12223 ] , which appeared at the time 
of submitting this paper. 
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